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O ■ ABSTRACT 

(~| I We build a simple, top-down model for the gas density and temperature profiles for 

'pHj galaxy clusters. The gas is assumed to be in hydrostatic equilibrium along with a com- 

O ' ponent of non-thermal pressure taken from simulations and the gas fraction approaches 

"^ ' the cosmic mean value only at the virial radius or beyond. The free parameters of 

^^' the model are the slope and normalisation of the concentration-mass relation, the gas 



polytropic index, and slope and normalisation of the mass-temperature relation. These 
parameters can be fixed from X-Ray and lensing observations. We coi npare our gas 



(N 

>, 
O ! pressu re profiles to the recently proposed 'Universal' pressure profile by lArnaud et al 

S^ I ( 2009 ) and find very good agreement. We find that the Sunyaev-Zel'dovich Effect (SZE) 

O ' scaling relations between the integrated SZE flux, Y, the cluster gas temperature, Tsi, 

10 ' the cluster mass, Mtot, and the gas mass, M^i. ^ are in excellent agree ment with the 

o ■ 

O 



recentl y observed r25n n SZE scaling relations bv lBonamente et al.l (J2008l l and rsoo rela- 



tion by lArnaud et al.l (120091 ). The gas mass fraction increases with cluster mass and is 



given by /gas(r50o) = 0.1324 + 0.0284 log ( iqis^-Tm. )• This is within 10% of observed 
k> i /gas (''500 )• The consistency between the global properties of clusters detected in X-Rays 

H ' and in SZE shows that we are looking at a common population of clusters as a whole, 

and there is no deficit of SZE flux relative to expectations from X-Ray scaling prop- 
erties. Thus, it makes it easier to compare and cross-calibrate clusters from upcoming 
X-Ray and SZE surveys. 

Subject headings: galaxies: clusters: general — cosmology: miscellaneous 



1. Introduction 

Large yield SZE cluster surveys promise to do precision cosmology once cluster ma ss-observable 



scaling relations are reliably calibrated. This can be done through cluster observations (JBenson et al 



2004: Bonamente et al 



modeling (JBulbul et al 



200S ) , simulations (|da Silva et alj l2004l : iBonaldi et alj l2007l ) and analytic 



20ld ). It is well known that different a strophysical process e s influence the 



cluste r mass-observable relati o ns non-trivially ( for ex ample, see lBalogh et alj (J2001h : lBorgani et al 



(|2004l ): lKravtsov et al.1 (J2005l ) : iPuchwein et al.1 (J2008l )) which can lead to biases in determining cos- 
mology with clusters. Alterii atively, one can 'self-calibrate' the uncertainties (jMaiumdar &: Mohr 



2003 



20041 : iLima fc Hull2004l V 



Simplistic modeling of the intra-cluster medium (ICM), like the 'isothermal /3- model' can give 
rise to inaccuracies. More complex modeling needs additional assumptions (such as gas following 



baryonic physics (jOstriker et al.ll2005l ) 



dark matter at la rge radii (IKoniatsu &: Seliakll200ll ). hereafter KS) or inclusion of less understood 



To partially circumvent our incomplete knowledge of cluster gas physics, we build a top-down 
phenomenological model of cluster structure, taking clues from both observations and simulations. 
It stands on three simple, well motivated, assumptions: (i) present X-Ray observations can give 
reliable cluster mass-temperature relations at r < r5oo which is used to cali brate our mode l s; (ii 



the gas rnass fr action, fgas, increase s with radius as s een in observations (jVikhlinin et al.l 12006 



Sun et al.ll2009l ) and in simulations (jEttori et al.ll2006l ). with non-gravitational processes pushing 
the gas outwards. It reaches values close to universal baryon fraction at or beyond the virial radius; 
and (iii) there is a component of non-thermal pressure support whose value relative to thermal 
pressure can be inferred from biases in mass estimates found in sim ulations (see 'Rasia et al.l ( 20041 )). 
This simple model can reproduce the 'Universal' pressure profile (lArnaud et al..,2009 ). X-Ray gas 
fraction, and SZE scaling relations in excellent agreement with observations. 



The Cluster Model 



2.1. The Cluster Mass Profiles 



The NEW profile (JNavarro et al.l 119971 ) is typically used to describe the dark matter mass 
profile. Here we adopt a NEW form for t he total matter profile since we use the observationally esti- 
mated concentration parameter given bv lComerford &: NataraianI (|2007l ) Cmr = ^'^f^^''^ (Myir/M^ 
Here M^ = 1.3 x 10^^ h~^MQ. The virial radius, rvir, is calculated from the spherical collapse model 

(|PeebleslE98ol l as r^ir = \ j^-^¥fir-r{] ■ Here, Ac{z) = ISvr^ + 82x - 39x2 (JBrvan fc Norman 



0.15±0.13 



-pcrit{z)Ac{z) 



19981 ) and x = fi„j(z) - 1. 



2.2. The Temperature and density Profiles 

XMM- Newton and Chandra observations have s hown that the cluster temperature declines at 



large radii (lArnaud et al. 



2005; 



Vikhlinin et al. 



20061) for both c ool (CC) and non-cool core (NCC) 



clusters. Simulations ( Ascasibar et al.ll2003l : iBorgani et al.ll2004l ). observations ( Sanderson &: Ponman 



20 id ) and analytic studies (JBulbul et al.1 I2OIOI ) indicate polytropic profiles for gas temperature. 
These studies also point towards an almost constant polytropic index 7 ~ 1.2 (atleast, till rsoo)- 
Hence, we adopt T{r) = T{0)f{r)'^-^ and p{r) = p(0)/(r) . We take the fiducial 7 = 1.2. 



We have also com pared the resulting temperature profiles with recent observations (jPratt et al 



20071 : ISun et al.ll2009l ) and find the decrements to be comparable . Alternatively, for our "Best-fit" 
models we let the 7 vary along with the Mq and a of the M-T relation and find the values that 
give the best-fit to the SZ Scaling relations. Further, CC clusters are charac terized by central 
temperature decrements which we take to be T{r) oc r' 



,0.3 



([Sanderson et al.l 120061 ) below O.lrsoo 



To calculate ICM d ensity and temperature profiles, we use the gas dynamical equation (JBinney &: Tremaine 



19871 : iRasiaelaLllaOOJ): 



dr 



1 dP{r) 1 



p{r) dr p{r) 



d[p{r)al{r)] a^ijr) 
- + 2p[r) 



dr 



(1) 



where $ is the gravitational potential, a'^{r) is the gas velocity dispersion, /3(r) is the velocity 
dispersion anisotropy parameter (put equal to zero in this work) and P{r) and p{r), the gas pressure 

P}LLh^(rr.\ Here, rup is the proton mass and // is the mean molecular 

1 dP{r) 



and density where P{r) = ^-^kT(r 



p{r) — «(«=- 



This is 



weight. For hydrostatic equilibrium without non-thermal pressure, 

normally used to obtain ICM profiles for a given halo (for example, in lKomatsu &: SeliakI (J200ll )). 



d$(r) 
I — efr' — 



Simulations (Rasia et al. 



2004 



Battaglia et al.ll2010l) show that non-thermal pressure ca n be 



significant es pecially at large radii. Both observ ations (JMahdavi et al.ll2008l : IZhang et al.ll2008l ) and 
simulations (JNagai et al.l 120071 : iLau et al.l l2009l ) suggest that the cluster mass calculated assuming 
only hydrostatic equilibrium is less than the true mass of a cluster. This discrepancy increases 
with radius. Typical values are 20-40% at rvir. The velocity dispersion term arises from the bulk 
motions of the ICM and contributes to the non-thermal pressure support. The profile f(r) can thus 
be numerically obtained by solving equation 1. 



2.2.1. Temperature and Density Normalization 
The temper ature profiles are normalized to the recently observed X-Ray M500 — T'sp scaling 



relation found by ISun et al.l (J2009l ) which includes data from cluster to group scales and is given by 



M^QoE{z) = Mo 



T. 



sp 



3keV 



(2) 



where Mq = (1.21 zb 0.08) x IO^^/i'^Mq and a = 1.68 ± .04. Here M500 is the mass within rsoo, 
where the average d ensity is SOOprJz ) whe re Pc[z) is the critical density at redshift z. Using the 
prescription given bv lMazzotta et al.l (J2004l ) we estimate the 'spectroscopic-like' temperature Tgi, a 
particular weighted average of T{r). This value of a 7^ 1.5 indicates deviation from self-similarity, 
pointing to non-gravitational energetics in the ICM. Here we bypass the microphysics that breaks 




0.01 - 



Fig. 1. — The normalized ICM pressure, P/Pmr is plotted against cluster radius for a cluster of 
mass 5 x 10^^ /i~^ Mq and z = 0. The thick black solid line is the fiducial model; the blue dot-dashed 
line is for Pnon-th = 0; the green dashed line has polytropic index changed to 1.12 from 1.2; the red 
dotted line is for lower concentration. The KS model pressure is given by the black solid line with 
circles. Note, that the SZE flux is given by the line-of-sight integral of P{r) over a given cluster 
area. 



'self-similarity' but normalize the cluster temperatures so as to exactly reproduce the observed 
M — Tsi relation. Thus our cluster model can be thought of as a top-down model. 

For any point in parameter space, representative of a simulated cluster, we calculate analyt- 
ically the temperature and density profiles. We start with an initial arbitrary T(0) and solve for 
f(r) as described earlier. Next, Tgi, is calculated in the radial range O.lrsoo - ?^5oo- The original 



T(0) is now adjusted by the ratio of Tgi to the Tgp from the observed M; 



^500 



Tgp relation. 



The 



equation for /(r) is now solved with this new T(0) after which the Tsi is again calculated. In a few 
iterations, a self consistent profile f(r) is obtained. Next, p{0) is determined by equating the fgas 
within the cluster radius to 0.9{Qb/^m) at the c luster boundary (r2nn or beyond). The Universal 
baryon fraction fift/Qm is given by 0.167 ± .009 (JKomatsu et al.ll20ld ). 



Simulations (JEttori et al.ll2006l ) and observations ()Vikhlinin et al.ll2006l : ISun et al.ll2009l ) show 



Mgas(r-) 



the gas mass fraction, fgas{'r) = m i-r) ' increases with radius. Stellar mass which accounts for a 
finite fraction of the baryons is larger at smaller radii such as ^2500 and for group scale haloes, as 
observed in the above mentioned studies. Radiative simulations tend to underestimate fgas due to 
overcooling and predict fgas = 0.7 — O.8(0b/r2m) at Vyir. Allowing 10% of the baryons to form stars, 
we take fgas = ^-^iS^b/^m) at the cluster boundary. The resulting fgas as seen in fig [2] shows good 
agreement with the observations at r5 nn. We assume that n on-gravitational effects o nly redistribute 



2010 



the gas. Recently, both observa tions (JRasheed et al.ll2010l ) and theoretical studies (JBattaglia et al. 



Nath &: Maiumdaijl2010l ) show that gas is driven outside the virial radius r^r and atleast 



10 , 10 

500 ^ sol' 




Fig. 2. — Left Panel: The fgas at rsoo is plotted against M500 for the cluster. The topmost 
horizontal line is for /gas = q^] the solid and dot-dashed lines are for KS model for two different 
halo concentr ations; the dotted lin e is for our fiducial model and the dashed line is from X-Ray 



observations (IVikhlinin et al.ll2009l ). Right P anel: The Y — Mp^s sca ling relations at r25oo for our 



models are shown with the data points from iBonamente et al.l ([2003) and the best fit Bonamente 
relation. The 5 lines pertaining to our models are - (i) thin blue solid line for model 1 (fiducial), 
(ii)thin blue dashed line for model 2 , (iii)thin blue dot-dashed line for model 3, (iv) Thick red 
dot-dashed line for the model 'Bestfit-1' and (v) thick red dashed line for the model 'Bestfit-3' (see 
text for details). The Bonamente bestfit to data is given by the black solid line with stars. The KS 
model is shown by dashed line with circles. 



upto 2r^ 



2.3. Model Descriptions 

We include non-thermal pressure, Pnon-th; in our calculations. However, this contribution to 
the total pressure {Ptot) for a cluster is difficult to model analytically. 



In this work, we follow gas dynamical simulations bv lRasia et al.l (J2004l ) to estimate the -Pnon-th- 
We adopt their Pth/Ptot as an input to our model. The mass of the cluster calculated from the 
hydrostatic term only is lower than the true mass by ~ 15% at rsoo, ~ 30% at r^ir and ~ 40% at 
2r^jr in our fiducial model. These values when compared with figure 13 in their paper are found to 
be of comparable magnitude. 

We consider the following models : 



model 1 (the fiducial model): Here /gas = 0.9(o-^) at r 



Mo = 1.728 X IOI^'Mq, a 




Fig. 3. — Left panel : The Y — T^i relation. The lines here have the same meaning as in figure [21 
except that the thick red dashed line is now for our model 'Bestfit-2' (see text). Right panel: The 
Y — Mtot relation at r25oo is plotted. The line scheme is the same. 



1.68 and 7 =1.2. We follow iRasia et al.l (|200J) to estimate -Pnon-th/^totai- 

• model 2: similar to model 1 but /gas = 0.9(fy^) at r = 2rvir. Pnon-th/-ftotai is extrapolated 
beyond r^r following simulations by Rasiqj- 

• model 3: parameters same as in model 1 but for 'zero' Pnon-th- 

Other than these models, we look at variations of the fiducial model, where we vary the parameters 
Mq, a and 7, to get the best fit to the lBonamente et al.l (J2008l ) SZE data. These are called Bestfit-1, 

2 2 2 



Bestfit-2 and Bestfit-3 and give a minimum to Xtot = xfy-T ) 
and xfy-M ) respectively, where the x^ is to Bonamente data. 

In figure [U we show the effect of varying some of the model parameters on the ICM pressure 
profile for a 5 x 10^^/i^^Mq cluster normali zed to P,„> which is taken to be the ICM pressure at 
r^ir for the standard self-similar model (see lArnaud et al.l (J2009l ) , Appendix A) . Inclusion of non- 



thermal pressure leads to shallower slope at large radii compared to only thermal pressure. The 
polytropic index has little infiuence on the pressure profile for the given change in 7. The integrated 
SZE, unlike X-Ray, is similar for both CC and NCC clusters. In figure [2] (left panel), we show that 
clusters in our model naturally have a mass dependent gas fraction, in agreement with observed 
/gas to within 10% for M^oq> 2 x 10^^h~^MQ. For our fiducial model, we find : 
/gas(r50o) = 0.1324 + 0.0284 log( ,oi,^^_°i° ). 



^Private communications 



Table 1: SZ Scaling Relations at r25oo : ^og{YD\E{z)^) = A + B\og{X/oj,) where ct,i = 
8keV, Cm^,, = 3* 10^3^0 and c^ntot = 3* IQi^M© and {5 = 1, -2/3, -2/3) for {X = T,i, Mg^s, Aftot) 
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Y-Mtot 
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Ax' 


A 
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Ax' 
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Ax' 


Bonamente 


-4.10±0.22 


2.37±0.23 


.0017 


-4.25±1.77 


1.41±0.13 


.073 


-4.20±3.00 


1.66±0.20 


.047 


Best fit 


-4.094 


2.363 




-4.25 


1.414 




-4.19 


1.654 




modell 


-4.215 


2.906 


0.57 


-4.258 


1.565 


0.12 


-4.322 


2.003 


0.71 


model2 


-4.342 


3.149 


2.04 


-4.201 


1.495 


0.024 


-4.441 


2.170 


1.83 


models 


-4.093 


2.944 


.19 


-4.376 


1.56 


1.31 


-4.245 


2.047 


0.24 


K-S 


-4.410 


2.28 


6.64 


-3.630 


2.18 


11.55 


-4.050 


1.95 


0.46 


Bestfit-1 


-4.153 


2.909 


0.20 


-4.301 


1.533 


0.39 


-4.247 


1.902 


0.23 


Bestfit-2 


-4.103 


2.448 


.0067 


- 


- 


- 


- 


- 


- 


Bestfit-3 


- 


- 


- 


-4.207 


1.544 


0.038 


- 


- 


- 



3. The SZE scaling relations - observations and and theoretical models 



The measurement of SZE (jSunyaev fc ZeldovichlllQSOl ) has come of age in recent times with im- 
provement in detector technologies. Both targeted observations (say from OVRO/ BIMA/SZA) and 



blank sky surveys (ACT/SPT) are underway having much cosmological potential (jCarlstrom et al. 



2002I ). Targeted observations have recently given us the SZE scaling relations which can now be 
used in surveys as proxy for mass. 



The SZE scaling relations (JBonamente et al.l 120081 ) predicted from self-similar theory are : 



YDi 



oc 



f T 



5/2 



E{z) 



YD\ oc /, 



M!l?E{zfl^ 



gas^^J-tot 



YD\ 



oc f^aTM^l^Ei 



n2/3 



gas 



(3) 



where Y is the the integrated SZE flux from the cluster, D^ is the angular diameter distance. Mg, 
and Mtot are the gas mass and total mass. 



3.1. The r25oo Scaling Relations 



Benson et al.l (J2004l ) presented the first observed SZE scaling relations between the c entral 



decrement, yo,Y and T^i for a sample of 14 clusters. Recently, iBonamente et al.l (J2008l ) have 
published scaling relations for 38 clusters at 0.14 < z < 0.89 using Chandra X-Ray observations 
and radio observations with BIMA / OVRO. Weak l ensing mass measurements, at ?'4ooo-8000) of 
SZE clusters have now been done by iMarrone et al.l (J2009l ) to give the Y — Mtot scaling. Their 
extrapolated masses at r25oo show agreement to within 20% to the hydrostatic mass estimates by 



Bonamente et alj (|2008l ). We follow I Bonamente et al.1 (|2008l ) in constructing our scaling relations. 
In particular we fit beta profiles to the density profile as well as the x-ray surface brightness Sx 
and compton y parameter obtained by projecting the temperature and density profiles obtained as 
a result of solving equation 1. : 

kBCTT 



nipC^ 



Tup dl: 



Sx oc 



Up dl 



(4) 



The isothermal temperature for each cluster is calculated in the same radial annulus as theirs. 
The SZE flux is found by integrating the SZE /3-profiles and the total mass assuming hydrostatic 
equilibrium is estimated using Mtot{r) = q (^^^IT^^n • With this prescription, we construct the 
three power law scaling relations given in equation [3l The coefficients for these scaling relations 
are specified in table [H For comparison, we a lso calculate the SZE sca ling relations from /3-fits to 
the often used 'Komatsu-Seljak' (KS) model (JKomatsu fc Seljaldl200ll ). A comparison of the SZE 



scaling relation for our models, the KS model and the Bonamente data are shown in figures [2] and 
[3] and table [TJ 

The first point to notice is the good agreement of our model scaling relations with the Bona- 
mente best fit. Especially, for the Y — Mgas relation, our models are in excellent agreement with 
observations. This is relevant as observationally the Y — Mgas relation has the least uncertainty. 
However the assumption of a /3- model for the ICM adds to the uncertainty. For the Y — Ts\ relation, 
our estimate of Tsi is not accurate for lower temperatures and the agreement with the Bonamente 
data becomes worse. Especially, for lower T^i, our models under predict the SZE flux. Note, that for 
both of these relations, the KS model line lies outside the data points and hence is a very bad fit to 
SZE scalings. Our models also under predicts the SZE flux for masses below M2500 ~ 3 x 10^^M(rc|. 

Next, we discuss our 'Bestfit' models. Once we vary the amplitude and slope of the M — Tgi 
relation and the polytropic index 7, the best fit values of these parameters obtained are in broad 
agreement with X-Ray observations. For example, for 'Bestfit-1', where we add the x^ from all the 
three scaling relations, our recovered values are (Mo,a) = (1.73 x 10^^Mq,1.7) which are within 
1 — o" of the X-Ray values (jSun et al.ll2009l ). The best fits are weakly sensitive to the value of 7; 
the 'Bestfit-1' model prefers a 7 = 1.14 which is lower than our fiducial value for 7. 

In general, the present data has large error bars and scatter and cannot distinguish between 
different models (with the exception of the KS model). However, there are three main points to 
note: (i) Y — Mgas is affected more by non-thermal pressure, since its presence influences how much 
gas can be pushed out. Our models are within 1 — o" of the Bonamente best fit for Y — Mgas while 
KS model is > 3 — o" away; (ii) At r2500j non-thermal pressure has lesser influence on the pressure 
support. Hence, the 'only thermal pressure' model is a good fit to the Y — Mtot data, followed 
by the fiducial model. Here, KS model, with no non-thermal pressure, is also within 1 — o" to the 
best fit; (iii) Since T^i is found by averaging over an region around ^25005 it is less infiuenced by the 



^The current limiting mass for both SPT and ACT surveys is higher. 



presence of non-thermal pressure and hence the trend in y — Tgi is similar to the trend in y — Mtot • 
However, KS model with its adiabatic normalization of T{r) is 2 — a away from Bonamente best 
fit. 



3.2. The r5oo Scaling Relations obtained from XRay Observations 

We compare our pressure profiles and scaling relatio r is wit h the recent 'Universal' pressure pro- 
file and resulting SZ scaling obtained by lArnaud et al.l (|2009l ) from X-Ray observat i ons. We also 
compare in fig. |3]the pressure profile with those obtained recently bv lBattaglia et al.l ( 201Cll). which 



comes from hydro simulations incorporating a prescription for AGN feedback, and ISehgal et al. 



( 2010l) where hot gas distribution within halos is calculated using a hydrostatic equilibrium model 
(JBode et al.l |2009| ). Between 0.1i?5oo — ^500 j i-e- the core radius and the upper limit for the 



X-Ray observations, all pressure profiles agree with the observations to within 20%. All the the- 
oretical pressure profiles start deviating significantly from the observed profile beyond ^500- From 
the pressure p rofile, we construct the scaling relation Y^qq = 10-^ (M500/3 x IO^'^/ijq ) hjQ . 
Arnaud et all (|20oi) find B = -4.739 ± 0.003 and A = 1.790 ± 0.015. We obtain {B,A) = 
(-4.646, 1.670) and (-4.797, 1.805) for Mode l 1 an d Mo del 2 respect i vely. For the sake of com- 
parison, the values found for iBattaglia et al.l (|2010l ) and ISehgal et al.l (|2010l ) pressure profiles are 
{B, A) = (-4.5 ± 0.1, 1.75 ± 0.06) and (-4.713 ± 0.004, 1.668 ± 0.009). 




Fig. 4. — Compari son of simulate d and semi-analytic gas pressure profiles with observed 'Universal' 
pressure profile by lArnaud et al.l (120091 ) for a cluster having M500 = 2 x lO^^/i~^M0. Plotted are 
the fractional differences of our fidu cial model given by red solid line, th e prof iles obtained from 
simulations bv lBattaglia et al.l (|20ld ) (gree n dashed line) a.ndlSehgal et al.l (|20ld ) (blue dot-dashed 
line) w.r.t the universal profile found by lArnaud et al.l (|2009l ) from observations upto rsoo and 
simulations beyond. 
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4. Discussions and Conclusion 

We have constructed a top-down model for galaxy clusters, normalized to the mass-temperature 
relation from X-Ray observations. The gas density and temperature profiles are found by iteratively 
solving the gas dynamical equation having both t hermal and non-th ermal pressure support. The 
form of the non-thermal pressure used is taken from lRasia et al.l ( 2004 ) . In our model, /gas becomes 
0.9 (q-^) at the cluster boundary , whereas gas is pushed out of the cluster cores to give /gas(?"50o) = 
0.1324 + 0.0284 log( ^Qiffio^^ ), similar to X-Ray observations. 



At r2500! the SZE scaling relations between SZE flux Y and the cluster average temperature, 
Tsi, gas mass, Mg as, and total mass , Mtpt , show very good agreement and are within 1 — cr to the 
best fit line to the lBonamente et al.l ( 20081 ) data. Especially, for the y — Mgas relation the agreement 



is excellent. In comparison, we also show that the Komatsu-Seljak model is in less agreement to 
the SZE scaling relations, especially for Y — T^i and Y — Mgas. Our ro ^ pn sca ling relations can be 
compared to those obtained from simulations. For example, the iNagai I (|2006l )(see their table 3 for 
scaling parameters) radiative simulati o n pre diction for the Y — Mgas relation giv es a Ax^ = 4.7 
w.r.t. to the best fit iBonamente et al.l ()2008l ) relation. Recently iBode et al.l (J2009l ) have predicted 
SZ scalings from a mixture of N-body simulations plus semi-analytic gas models, normalized to 
X-Ray observations for low-z clusters. The Ax^ of their model is 0.06 and agrees well with our 
results. 

Further out, at t^qq, the Y — M scali ng relation obtained for our models agree very well with 
those obtained from X-Ray observations (jArnaud et al.ll2009l ). Most assuringly, the gas pressure 
profile in our simple phenomenological model of clusters, comes out to be with in ~ 20% beyond .1 
R500 to the observed 'Universal' pressure profile given by lArnaud et al.l (J2009l ). 



Most importantly, the fact that X-Ray normalised models can reproduce SZE scaling relations 
well is reassuring for cluster studies. It shows that we are looking at a common population of 
clusters as a whole, and there is no deficit of SZE flux relative to expectations from X-Ray scaling 
properties. Thus, one can compare and cross-calibrate clusters from upcoming X-Ray and SZE 
surveys with increased confidence. It also gives us confidence to extrapolate our models to larger 
radii in order to construct the Y — M200 scaling relation and SZ power spectrum templates. 
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